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We  present  the  application  of  dynamic  programming  to  a  combinatorial  optimization 
problem  to  achieve  proper  arrival  runway  spacing,  which  appears  in  the  process  of  assign¬ 
ing  speed  during  the  transition  to  approach  and  approach  phases  of  flight.  We  apply  the 
algorithm  to  data  from  a  fast-time  simulation  developed  under  NASA’s  Advanced  Air 
Transportation  Technologies  Project  for  investigating  new  air  traffic  management  (ATM) 
concepts.  For  this  research,  the  simulation  is  configured  to  simulate  traffic  inbound  along 
two  arrivals  to  the  Dallas/Ft.  Worth  (DFW)  airport,  merging  into  a  single  stream  at  fix 
just  prior  to  the  final  approach  fix.  We  show  how  the  algorithm  computes  the  maximum 
minimum  spacing  between  aircraft  upon  landing,  and  investigate  the  sensitivity  of  the 
spacing  to  perturbations. 
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Introduction 

The  recent  growth  in  embedded  systems  control 
applications  has  been  triggered  by  numerous  factors, 
which  range  from  advances  in  monitoring  technology, 
communications,  and  available  computing  power.  It 
has  provided  engineers  unprecedented  capabilities  for 
real-time  control.  Yet  a  bottleneck  still  exists  in  the 
design  and  realization  of  control  laws  for  large  scale 
systems,  due  to  the  complexity  of  the  underlying  math¬ 
ematical  problems.  In  particular,  NP-completeness 
or  non  convexity  of  optimization  problems  cannot  be 
eased  by  increases  in  computational  power  alone,  but 
by  new  breakthroughs  in  algorithm  design.  Indeed, 
this  type  of  difficulty  does  not  preclude  the  tractabil- 
ity  of  solutions  which  are  close  to  the  optimal,  fast  to 
compute,  and  guaranteed  to  achieve  a  certain  perfor- 
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mance.  Combinatorial  optimization  algorithms  have 
proved  to  be  very  useful  to  solve  this  type  of  difficul¬ 
ties,  and  have  emerged  as  a  very  powerful  tool  to  solve 
planning  problems  for  large  scale  systems,  such  as  Air 
Traffic  Control  (ATC),  or  networks  of  multiple  vehi¬ 
cles. 

Several  recent  publications  have  used  efficient  al¬ 
gorithmic  techniques  to  solve  combinatorial  optimiza¬ 
tion  problems  which  appear  naturally  in  large  scale 
systems.  In  particular,  for  ATC,  Ribichini  and  Fraz- 
zoli  [10]  reduce  aircraft  coordination  problems  to  an 
integer  program  similar  to  facility  location,  for  which 
they  develop  a  decentralized  approximation  algorithm. 
Neogi  [7]  reduces  the  problem  of  landing  aircraft 
assignment  for  multiple  runway  configurations  to  a 
Steiner  Tree  problem.  Parallel  landing  with  aircraft 
dependent  separation  requirement  is  known  to  be  an 
NP-complete  scheduling  problem  [8].  In  [4],  we  de¬ 
velop  a  polynomial  time  algorithm  to  solve  the  prob¬ 
lem  of  maximum  spacing  between  aircraft  when  the 
set  of  feasible  times  of  arrival  is  a  single  interval  for 
every  aircraft,  and  the  landing  order  is  not  known  a 
priori.  We  also  developed  an  approximation  algorithm 
for  minimizing  the  latest  landing  time  of  a  platoon  of 
aircraft,  the  case  in  which  the  aircraft  can  hold  for 
periodic  time  intervals  [2].  This  list  of  problems  is 
not  exhaustive,  and  becomes  very  large  if  Unmanned 
Aerial  Vehicle  (UAV)  problems  are  included.  In  this 
paper,  we  will  be  interested  in  ATC  applications. 

The  ATM  system  is  a  large-scale  complex  system 
that  has  evolved  over  decades  to  become  highly  ro- 
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bust  and  safe.  It  includes  both  human  agents,  such 
as  pilots,  air  traffic  service  providers  (ATSPs),  and 
airline  operations  center  (AOC)  personnel,  and  ma¬ 
chine  agents  (e.g.,  aircraft  flight  management  systems 
(FMSs)  and  ATC  automation).  The  present-day  ATM 
system  remains  vulnerable  to  environmental  distur¬ 
bances  and  unpredictable  in  many  respects.  It  relies 
on  conservative  margins  and  closely  monitored  tacti¬ 
cal  operations  to  ensure  safety,  which  in  turn  leads  to 
inefficiencies. 

Researchers  seek  a  future  ATM  system  that  lever¬ 
ages  advances  in  aircraft  automation,  ATC  automa¬ 
tion,  and  communications,  navigation,  and  surveil¬ 
lance  systems  to  decrease  delays,  noise,  and  fuel  usage 
while  improving  safety.  One  class  of  proposed  future 
ATM  concepts  entails  keeping  aircraft  on  FMS  routes 
until  stabilized  on  final  approach,  enabling  them  to  fly 
more  precise  and  efficient  trajectories  in  the  terminal 
(TRACON)  airspace.  For  example,  FMS  procedures 
enable  ’continuous  descent  approaches  (CD As),’  which 
afford  greater  fuel  efficiency  and  reduced  noise  if  they 
can  be  flown  safely. 

Efficient  TRACON  ATM  rests  on  the  ATSP’s  ability 
to  predict  future  aircraft  locations  and  issue  clearances 
appropriate  for  merging  aircraft  and  maintaining  re¬ 
quired  separation  and  flow  spacing.  TRACON  ATSPs 
use  speed  changes  and  lateral  maneuvers  to  control 
traffic.  Control  options  are  limited  by  both  the  area 
available  for  maneuvering  and  the  overriding  concern 
of  pilots  to  stabilize  the  aircraft  for  the  approach.  Ac¬ 
celerations  are  therefore  not  likely  to  be  desirable  to 
pilots,  while  lateral  maneuvers  may  lead  to  adjusting 
the  overall  arrival  flow  if  traffic  density  is  sufficiently 
high. 

Studies  in  today’s  ATM  environment  have  shown 
that  ATSPs  find  it  difficult  to  predict  the  4D  trajec¬ 
tory  of  aircraft  flying  CDAs.  ATSPs  therefore  add  a 
large  buffer  to  ensure  safety,  which  dramatically  limits 
throughput  [7].  One  challenge  is  therefore  to  develop 
an  ATM  concept  that  yields  CDA-type  benefits  with 
moderate  to  high  throughput.  ATC  automation  could 
potentially  support  ATSPs  toward  this  end.  For  exam¬ 
ple,  automation  tools  could  use  4D  FMS  trajectories 
downlinked  by  aircraft  to  provide  ATSPs  with  advi¬ 
sories  on  suitable  speed  changes  to  issue  aircraft  to  null 
schedule  errors  created  by  environmental  factors.  In 
this  paper,  we  address  issues  surrounding  speed  advi¬ 
sories  used  in  conjunction  with  scheduling  automation 
to  reduce  excess  arrival  spacing  buffers. 

Investigating  new  ATM  concepts  is  inherently  dif¬ 
ficult  due  to  the  number  of  agents  and  overall  com¬ 
plexity  of  the  system.  Human  agents  are  essential 
for  system  robustness,  and  acceptability  of  procedures 
and  automation  tools  to  human  agents  is  crucial.  A 
promising  concept  must  therefore  undergo  extensive 
human-in-the-loop  simulation  to  refine  these  elements. 
However,  a  number  of  human  factors  issues  may  be  ex¬ 


Fig.  1  Arrival  chart  showing  merging  FMS  routes. 
Aircraft  are  expected  to  fly  the  arrival  in  LNAV/VNAV 
nntil  glideslope  intercept. 


amined  using  simulations  with  embedded  agent  mod¬ 
els.  These  include  the  capability  of  ATC  automation 
to  generate  usable  advisories,  the  pace  at  which  ATSPs 
must  issue  various  types  of  clearances  to  control  traffic 
under  various  conditions,  the  types  and  magnitudes  of 
disturbances  that  ATSPs  can  manage  given  proposed 
routings  and  clearances,  and  related  issues.  In  this 
research  we  use  a  Java^'^'^-based  fast-time  simulation 
tool  called  TCSim  (Trajectory-Centered  Simulation) 
to  provide  relevant  data  and  visualize  proposed  oper¬ 
ational  concepts. 

This  paper  is  organized  as  follows.  We  first  describe 
the  capabilities  of  TCSim  and  how  we  use  it  to  sim¬ 
ulate  a  specific  arrival  problem  at  the  DFW  airport. 
Then,  we  describe  combinatorial  optimization  prob¬ 
lems  which  appear  in  TCSim,  and  derive  a  polynomial 
time  algorithm  to  solve  it,  using  dynamic  program¬ 
ming.  We  finally  show  results  of  simulations  run  for 
the  DFW  airport,  which  illustrate  the  use  of  TCSim 
to  study  the  sensitivity  of  spacing  to  disturbances. 

Spacing  arrival  traffic 

TCSim  was  developed  at  NASA  to  provide  capabili¬ 
ties  useful  for  investigating  ’trajectory-oriented’  ATM 
concepts.  We  are  using  it  here  to  investigate  a  prob¬ 
lem  that  involves  merging  arrival  flows  in  the  DFW 
TRACON  (see  Figure  1).  In  this  configuration,  north¬ 
west  and  southwest  arrivals  enter  the  TRACON  over 
the  BAMBE  and  FEVER  meter  fixes,  respectively. 
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merge  at  GIBBI,  and  land  on  runway  18R  spaced  ac¬ 
cording  to  a  standard  wake  vortex  spacing  matrix. 
This  problem  involves  several  issues,  including  suitable 
FMS  routes  that  aircraft  can  fly  using  typical  Lat¬ 
eral/Vertical  Navigation  (LNAV/VNAV)  modes,  air¬ 
craft  scheduling  at  the  meter  Axes  and  runway,  and 
the  design  of  control  strategies  and  supporting  ground 
automation  that  can  null  schedule  errors  to  produce 
an  efficient  arrival  flow.  Variations  and  extensions 
for  future  study  include  the  simultaneous  use  of  run¬ 
way  13R,  aircraft  equipped  to  meet  required-times-of- 
arrival  (RTAs)  assigned  by  ATSPs,  and  relative  self¬ 
spacing  guidance. 

TCSim  is  designed  to  simulate  common  types  of 
commercial  aircraft  flying  FMS  trajectories  in  fast 
time.  It  can  use  the  trajectories  it  generates  for 
each  aircraft  to  emulate  complex  air-  and  ground  au¬ 
tomation  functionality.  TCSim  also  includes  embed¬ 
ded  ATSP  and  pilot  agents  for  simulating  interac¬ 
tions  involved  with  specific  concepts  and  traffic  con¬ 
ditions.  TCSim  is  capable  of  automatically  generating 
traffic  scenarios  with  specific  characteristics  to  sup¬ 
port  Monte  Carlo  studies,  and  of  producing  a  vari¬ 
ety  of  metrics  for  each  trial.  Researchers  can  specify 
charted  routes,  traffic  characteristics,  agent  behavior, 
and  structure  of  simulation  trials  to  conduct  via  a  con¬ 
figuration  file. 

TCSim  represents  trajectories  as  a  sequence  of  lat¬ 
eral  legs  (including  turn  segments),  and  associated 
’vertical’  legs  that  specify  when  computed  altitude  and 
speed  rates  of  change  are  in  force.  A  rule  set  guides 
how  TCSim  constructs  vertical  trajectory  segments. 
TCSim  produces  vertical  trajectories  that  map  closely 
to  observed  VNAV  profiles.  Figure  2  shows  TCSim 
vertical  segment  types,  and  how  they  might  appear  in 
an  aircraft’s  approach  and  landing  trajectory.  Work¬ 
ing  backwards  from  touchdown,  each  aircraft  has  a 
segment  that  represents  a  stable  landing  configuration 
at  a  landing  speed  appropriate  for  the  aircraft’s  type. 
Prior  to  stabilization,  other  legs  along  the  glideslope 
represent  decelerations  for  successive  flap  extensions. 
The  glideslope  intercept  point  depends  on  the  crossing 
restriction  at  the  final  approach  fix;  aircraft  with  high 
landing  speeds  may  extend  flaps  earlier  than  shown  in 
Figure  2.  TCSim  constructs  trajectory  segments  that 
honor  the  type  (i.e.,  at,  at-or-above,  or  at-or-below) 
and  values  for  speed  and  altitude  restrictions.  Decel¬ 
erations  to  meet  crossing  speeds  use  a  deceleration  seg¬ 
ment  with  a  two  degree  flight  path  angle  in  the  absence 
of  other  constraints;  TCSim  can  also  model  steeper 
speedbrake-assisted  decelerations.  Crossing  restriction 
locations  dictate  the  need  for  trajectory  segments  with 
computed  flight  path  angles.  It  is  important  that 
charted  speed  and  altitude  restrictions  are  designed 
so  that  VNAV  can  reasonably  fly  them  in  the  given 
wind  held;  otherwise  TCSim  will  construct  unrealis¬ 
tic  trajectories.  Figure  3  plots  the  actual  trajectories 


TCSim  produces  for  aircraft  flying  the  charted  FMS 
arrival  routes  shown  in  Figure  1. 

Once  TCSim  computes  a  trajectory  for  a  particu¬ 
lar  aircraft,  it  ’flies’  the  trajectory  in  the  background 
to  collect  accurate  estimated  times-of-arrival  (ETAs) 
at  each  lateral  point  along  the  trajectory.  These 
estimates  support  emulation  of  automatically  down¬ 
linked  trajectory  information  for  input  to  ground- 
based  scheduling  automation.  TCSim  can  also  eval¬ 
uate  potential  changes  to  the  trajectory  and  supply 
ETAs  for  aircraft  given  the  changes.  This  capability 
is  essential  to  the  present  research;  we  use  TCSim  to 
output  a  set  of  feasible  arrival  times  for  each  arriving 
aircraft  that  serve  as  input  to  the  dynamic  program¬ 
ming  algorithm. 

TCSim  also  contains  embedded  ATSP  agents.  Sim¬ 
ulated  ATSPs  can  issue  a  variety  of  clearances  to 
aircraft,  including  heading,  altitude  and  speed  clear¬ 
ances,  as  well  as  more  complex  clearances  with  crossing 
restrictions,  and  a  variety  of  execution  conditions.  Re¬ 
gardless  of  clearance  type,  the  aircraft’s  flight  path  is 
still  represented  as  a  trajectory.  Pilot  agents  are  cur¬ 
rently  modeled  by  the  time  it  takes  and  the  method 
used  to  implement  a  particular  clearance  (TCSim  may 
also  simulate  flight  technical  errors,  including  clear¬ 
ance  non-compliance).  For  example,  the  new  trajec¬ 
tory  for  an  aircraft  cleared  to  fly  a  particular  course 
contains  a  short  delay  leg  at  the  beginning  to  represent 
the  aircraft’s  trajectory  during  the  time  the  crew  takes 
to  begin  the  turn.  After  the  delay  leg  the  trajectory 
includes  a  turn  segment  and.  Anally,  a  straight  seg¬ 
ment  on  the  specified  course.  Speed  clearances  operate 
similarly;  they  may  include  a  delay  segment  that  repre¬ 
sents  the  time  it  takes  a  crew  to  configure  the  aircraft 
automation  or  flight  surfaces.  Speed  clearances  that 
effectively  cancel  one  or  more  speed  restrictions  entail 
backwards  computation  of  the  trajectory  without  the 
cancelled  restrictions,  as  well  as  forward  integration 
through  the  wind  held  to  construct  the  required  speed 
change  segments. 

Each  ATSP  agent  uses  a  particular  ’control  strategy’ 
that  specifies  what  ATC  automation  is  available  to 
them  and  when  they  address  each  aircraft.  No  attempt 
is  made  to  simulate  the  actual  process  of  controlling 
air  traffic,  as  in  other  research  [5,  9].  Instead,  each 
controller  addresses  an  arriving  aircraft  as  it  passes 
specified  ’control  points’  (in  general,  the  farther  a 
control  point  is  from  the  runway,  the  more  potential 
control  authority  it  offers).  The  ATSP  agents  con¬ 
sult  simulated  ATC  automation  to  produce  specified 
advisories  for  that  aircraft.  When  actively  control¬ 
ling  traffic,  agents  can  select  an  advisory  to  null  each 
aircraft’s  arrival  schedule  error,  and  issue  it  as  a  clear¬ 
ance. 

TCSim  ATSP  agents  have  been  used,  for  exam¬ 
ple,  to  show  that  a  single  speed  advisory  clearance 
can  null  sufficiently  small  schedule  errors  caused  by 
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errors  in  meter  fix  crossing  (TRACON  entry)  times 
and  predicted  landing  speeds  [6].  In  that  research, 
the  TRACON  merge  problem  was  investigated  Monte 
Carlo-style  as  follows.  TCSim  automatically  gener¬ 
ated  arrival  flows  across  both  meter  fixes  to  runway 
18R.  It  first  generated  a  specified  number  of  aircraft  of 
types  selected  to  represent  a  typical  DFW  arrival  flow. 
It  assigned  each  aircraft  to  arrive  from  a  given  direc¬ 
tion  (based  on  a  specified  minimum  number  that  must 
come  from  each  direction),  and  generated  a  trajectory 
corresponding  to  the  appropriate  FMS  route.  TCSim 
then  established  simulation  entry  times  for  each  air¬ 
craft  that  simulated  a  coordinated  schedule  for  flows 
across  both  meter  fixes  to  yield  proper  wake  vortex 
spacing  at  the  runway  threshold.  It  then  adjusted 
the  entry  times  probabilistically  to  produce  meter  fix 
crossing  time  errors,  and  adjusted  the  landing  speeds 
and  recomputed  the  trajectory  to  simulate  predicted 
landing  speed  errors. 

With  the  scenario  specified,  TCSim  began  simu¬ 
lating  the  traffic  in  real-time.  ATSP  agents  issued 
automatically  generated  advisories  at  specified  control 
points  in  an  attempt  to  null  any  schedule  deviation 
for  each  aircraft.  We  specified  that  only  decelera¬ 
tions  could  be  issued,  on  the  assumption  that  pilots 
are  not  likely  to  find  accelerations  along  the  transition 
to  final  approach  desirable.  On  successive  trials  with 
each  traffic  scenario,  TCSim  incremented  an  ’excess 
spacing  buffer’  that  was  added  during  the  scheduling 
process.  When  the  buffer  was  small,  separation  vi¬ 
olations  at  the  merge  point  (GIBBI)  and  subsequent 
locations  were  observed.  The  results  of  multiple  tri¬ 
als  yielded  (among  other  useful  measures)  the  excess 
spacing  buffer  required  to  ensure  safe  operations  for 
the  given  traffic  scenario  under  the  given  control  strat¬ 
egy  with  disturbances  of  particular  types  sampled  from 
specific  distributions.  Results  of  this  sort  not  only 
indicate  the  effectiveness  of  a  given  control  strategy, 
but  also  inform  the  development  of  automation  tools, 
traffic  scenarios,  clearance  phraseology,  and  other  el¬ 


ements  required  to  conduct  useful  human-in-the-loop 
studies  of  future  ATM  concepts. 

In  the  research  reported  here,  we  use  TCSim  to  pro¬ 
vide  a  set  of  feasible  times  of  arrival  for  each  aircraft. 
The  times  are  represented  in  seconds  from  a  global 
reference  (elapsed  seconds  from  simulation  start).  The 
arrival  times  reflect  when  each  aircraft  will  arrive  given 
its  current  trajectory,  as  well  as  when  it  would  arrive 
if  issued  a  range  of  speed  clearances  in  five  knot  decre¬ 
ments  down  to  the  final  approach  fix  crossing  speed. 
The  algorithm  presented  below  can  use  this  informa¬ 
tion  to  automatically  space  arrival  traffic.  The  set  of 
feasible  arrival  times  provides  the  possibility  of  opti¬ 
mizing  the  spacing  according  to  a  user-defined  cost 
(which  is  arbitrary)  while  maintaining  safety.  For  the 
present  paper,  we  will  maximize  the  minimum  time 
separation  between  successive  landings. 

A  polynomial  time  algorithm 
for  maximal  spacing 

The  problem  of  computing  the  maximal  minimum 
spacing  of  jobs  scheduled  in  a  single  processor,  which 
we  investigate  in  this  paper,  is  not  a  standard  prob¬ 
lem  in  scheduling.  Usual  scheduling  problem  in  the 
literature  try  to  minimize  a  cost  relevant  for  computer 
science  or  operations  research  applications,  for  exam¬ 
ple  the  sum  of  all  processing  times,  the  makespan  (time 
at  which  the  last  job  is  finished),  and  the  sum  of  the 
delays  [1].  In  the  present  context,  we  are  interested 
in  maximizing  the  separation  between  landing  times, 
given  hard  constraints  which  define  the  set  of  possible 
arrival  times.  This  problem  was  posed  in  its  general 
form  as  a  mixed  integer  linear  program  (MILP)  in  [3], 
and  was  solved  for  specific  cases  in  [4,  2].  In  this  ar¬ 
ticle,  we  investigate  a  subcase  of  the  general  problem 
which  appeared  while  using  TCSim,  which  is  a  purely 
discrete  version  of  the  general  MILP  ^  in  [3] . 

^Note  that  in  general,  nothing  enables  to  know  a  priori  if 
the  transformation  from  a  MILP  into  an  integer  program  (IP) 
makes  it  easier  or  harder  to  solve. 
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Fig.  3  TCSim  altitude  and  speed  profiles  for  aircraft  flying  the  BAMBE  and  FEVER  FMS  arrivals  to  DEW 
runway  18R. 


Problem  formulation 

Consider  N  aircraft  labelled  by  the  index 
-  ,N}.  The  arrival  time  of  aircraft  i 
belongs  to  the  a  set  of  Ui  feasible  arrival  times 
(computed  by  TCSim),  called  tij,  indexed  by 
jS  {1,  •  •  •  ,ni}.  In  most  of  the  cases  run  in  TCSim, 
the  following  condition  holds:  for  all  i  and  t'  such  that 
i  >  i' ,  tij  >  ti'j!  for  all  j  and  j' .  This  means  that 
no  matter  what  choice  the  algorithm  makes,  aircraft 
always  arrive  in  increasing  order  of  i.  This  is  the  case 
whenever  the  perturbations  from  the  schedule  TCSim 
uses  to  generate  the  traffic  scenario  are  sufficiently 
small,  and  the  control  points  and  speed  restrictions 
along  the  arrival  route  are  selected  so  the  speeds  fall 
in  a  reasonably  narrow  range.  (The  TCSim  scenarios 
used  here  reflect  these  criteria:  perturbations  do 
not  exceed  approximately  thirty  seconds,  and  the 
available  speed  range  is  roughly  forty  knots.)  While 
ATSPs  could  conceivably  issue  clearances  (e.g., 
heading  vectors)  that  effectively  resequence  the  ar¬ 
rivals,  the  simulated  arrival  concept  requires  suitably 
small  deviations,  because  it  specifically  examines  the 
potential  for  using  speed  clearances  alone  to  properly 
space  the  aircraft.  For  this  paper,  we  want  to  find 
the  maximum  minimum  spacing  between  the  aircraft 
at  the  destination  airport,  mathematically  defined  by 
the  following  optimization  program: 

max:  A 

s.t.:  n  G  VjG  {1,  •  •  •  ,  iV}  (1) 

n  -  ri_i  >  A  ViG  {2,  •  •  •  ,  N} 

In  the  previous  program,  A  represents  the  maximum 
minimum  spacing  between  aircraft  upon  landing;  the 
ti,j  G  R  are  given  (computed  by  TCSim).  This  pro¬ 
gram  is  a  fully  discrete  version  of  the  MILP  solved 
in  polynomial  time  in  [4].  The  optimization  program 
solved  in  [4]  is  obtained  from  (1)  by  allowing  any  ar¬ 
rival  time  between  the  earliest  and  the  latest: 

rii 

Ti  G  [min  t ^ ,  max  ,  ] 
i=i  ’ 


and  allowing  any  order  of  arrival  of  the  aircraft,  i.e.: 

\n  -  Ti-i\  >  A 

Dynamic  programming  algorithm 

We  now  show  how  to  solve  this  problem  exactly  us¬ 
ing  dynamic  programming. 

Proposition  1.  The  solution  A  to  the  optimiza¬ 
tion  program  (1)  is  given  hy  A  =  (5(A,  uat),  where, 
5{N,nN)  is  defined  by  the  following  recurrence: 

<5(2,  j)  =  t2j  —  ti,i 

(5(i,  j)  =  niax{min{tjj  -  ,  i5(i-l,/)}} 

f  =1 

(2) 

The  computational  complexity  of  the  algorithm  is 

N 

0{TM)  where  T  =  is  the  total  number  of  ar- 

N 

rival  time  points,  and  M  =  maxTi^. 

Proof.  For  all  jG  {1,  •  •  •  ,  Ui},  <5(2,  j)  is  by  construction 
the  largest  separation  between  aircraft  1  and  2  if  air¬ 
craft  2  arrives  at  t2j,  since  tip  is  the  earliest  possible 
arrival  time  of  aircraft  1. 

Suppose  that  S{i—  1,/)  represents  the  largest  min¬ 
imum  separation  between  consecutive  aircraft  up  to 
i  —  1,  if  aircraft  i  —  1  arrives  at  time  ti-iji.  Then, 
min  {tij  —  ti-ij^  ,  S{i  —  l,j')}  is  the  minimum  be¬ 
tween  two  quantities:  (t)  the  separation  fip  —  tj-ij' 
between  aircraft  i  and  aircraft  i  —  1  if  aircraft  i  is 
scheduled  at  Uj  and  aircraft  i  —  1  at  ti-ip>;  and  (u) 
the  minimum  separation  S{i  —  l,j')  between  consec¬ 
utive  aircraft  up  to  j  —  1,  if  aircraft  i  —  1  arrives  at 
ti-ipi.  The  minimum  of  (t)  and  (tt)  is  the  minimum 
separation  between  consecutive  aircraft  up  to  i.  Tak¬ 
ing  the  maximum  over  the  set  {1,  •  •  •  ,ni}  provides  the 
largest  minimum  separation  between  the  i  consecutive 
aircraft. 

The  best  schedule  is  obtained  if  the  last  aircraft  ar¬ 
rives  as  late  as  possible,  i.e.  if  tn  =  tM,nN-  This 
number  is  6{N,n]s[) 
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There  are  T  evaluations  of  (5(-,  •)  to  do.  Each  of  them 
requires  taking  the  maximum  between  M  numbers. 
Therefore,  the  complexity  is  0{TM).  □ 

Corollary  2.  The  previous  result  holds  even  when  the 
arrival  times  overlap,  provided  the  recurrence  is  re¬ 
placed  by: 

Hhj)  = 

max  {minjtij  - 

/  G  {I,--  -  ,n*_i} 
ti  —  ^ 

Proof.  The  proof  is  the  same  as  for  Proposition  1.  For 
a  given  choice  of  Uj  for  aircraft  i,  the  set 

{j  \  j  G  {1,  *  *  *  ,  ,  ti—i  j'  ^  tij } 

represents  the  set  of  available  arrival  times  of  aircraft 
i  —  l  which  precede  ti  j,  and  therefore  respect  the  order 
of  arrival  of  the  aircraft.  □ 

Remark.  The  computation  6{i,j)  of  (2)  can  be  done 
in  log2(M)  time  using  a  bisection  method  on  j'  =  1  to 
Ui-i-  This  is  because  tij  —  ti-iji  is  decreasing  in  j' 
and  6{i—  1,  j')  is  increasing  in  f ,  so  that  the  maximum 
is  achieved  when  the  two  are  crossed  over.  This  can  be 
done  using  bisection  on  Ui-i  possible  points.  There¬ 
fore,  the  complexity  of  the  algorithm  is  0(riog2(M)). 
Note  that  M  is  small  in  practice;  therefore,  this  obser¬ 
vation  will  not  be  used  in  the  implementation. 

The  present  algorithm  requires  arrival  times  to  be 
known  a  priori.  The  problem  of  unknown  arrival  order 
might  conceivably  be  addressed  through  a  combination 
of  trajectory-based  prediction  and  downlinked  wind 
updates,  so  that  earlier  arrivals  enable  the  algorithm 
to  be  applied  to  later  arrivals.  Applications  of  this  sort 
are  beyond  the  scope  of  this  paper.  We  believe  that  a 
generalization  of  Baptiste’s  algorithm  [1]  would  solve 
this  problem  in  polynomial  time  (i.e.  even  if  the  order 
of  aircraft  is  not  fixed  anymore). 

Application  to  fast-time  simulation  data 

Application  to  fast-time  simulation  data.  The  algo¬ 
rithm  was  applied  to  TCSim  data  for  the  merging  flows 
at  DEW.  We  used  it  for  two  purposes.  First,  in  order 
to  find  the  maximum  minimum  spacing  upon  landing 
under  uncertainty.  Second,  in  order  to  show  how  to 
use  this  tool  to  quantify  the  perturbation  threshold 
above  which  it  is  not  possible  anymore  to  guarantee  a 
given  minimum  spacing  upon  landing.  Typically,  the 
type  of  data  generated  by  TCSim  reads  as  follows: 

ATAOOl  1295  1305  1310  1311  1316  1320  1325  1327  1335 
UAL002  1413  1423  1429  1439  1447  1458  1468  1478  1489 
DAL003  1522  1532  1541  1551  1557  1567  1581  1587  1601 
UAL004  1606  1613  1619  1629  1638  1648  1659  1673  1682 
COA005  1693  1700  1705  1710  1710  1715  1720  1725  1730 
SWA006  1787  1794  1799  1799  1804  1809  1814  1819  1824 
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where  the  names  in  the  first  column  are  dummy,  and 
the  arrival  time  is  in  seconds,  measured  from  an  ar¬ 
bitrary  reference.  We  realized  a  set  of  simulations  in¬ 
volving  20  aircraft  of  the  same  type.  This  assumption 
enables  us  to  ignore  differences  in  required  wake- vortex 
spacing  distances  due  to  different  aircraft  types. 

To  produce  the  data,  TCSim  performs  the  process 
described  above.  It  formulates  each  arrival  traffic  sce¬ 
nario  based  on  a  schedule,  but  produces  estimated 
arrival  times  that  reflect  TRACON-entry  time  and 
predicted  landing  speed  perturbations.  As  aircraft 
arriving  from  the  north  cross  HIKAY,  and  southern 
arrivals  cross  DELMO,  TCSim  computes  the  set  of 
feasible  arrival  times.  For  the  application  concept 
(i.e.,  speed  adjustments  alone),  we  assume  that  by  this 
point  the  ATSPs  will  have  assigned  the  landing  order 
for  the  merged  flow.  However,  the  exact  landing  time 
has  not  been  assigned. 

Maximum  minimum  spacing  upon  landing 

In  order  to  illustrate  the  variation  of  the  maximum 
minimum  spacing  with  the  number  of  aircraft  which 
have  already  crossed  this  point  (at  which  ATC  would 
issue  the  corresponding  speed  assignment),  we  com¬ 
pute  its  average  value  over  20  runs.  In  the  final  version 
of  the  paper,  we  will  display  the  result  of  a  Monte- 
Carlo  simulation.  Figure  5  shows  the  variation  of  A 
with  the  number  n  of  aircraft  which  have  crossed  the 
metering  point.  As  expected,  this  minimum  separation 
decreases  with  the  number  of  aircraft  (the  more  air¬ 
craft  are  taken  by  the  algorithm,  the  more  the  margin 
for  adjustment  decreases).  The  asymptote  is  around 
90  seconds,  which  is  around  the  value  expected  for  this 
type  of  data  set. 

Figure  6  illustrates  six  scenarios  for  which  we  have 
compared  the  maximum  minimum  spacing  provided 
by  the  algorithm  and  the  spacing  provided  by  TCSim. 
By  construction,  the  algorithm  does  not  take  other 
factors  than  time  spacing  into  account,  which  explains 
the  differences  between  the  two  results,  illustrated  in 
Figure  6.  Several  interesting  results  appear  from  these 
simulations: 

•  An  appearance  of  oscillations  is  observed  in  the 
TCSim  data.  This  is  in  fact  due  to  the  random 
nature  of  the  perturbations  from  TCSim  schedule, 
coupled  with  the  fact  that  the  simulated  ATSP 
control  strategy  only  allows  decelerations.  More¬ 
over,  for  this  research,  the  ATSP  agents  control 
every  aircraft  to  the  original  schedule,  including 
the  first  aircraft.  This  means  that  if  the  first 
aircraft  happens  to  be  ahead  of  schedule,  it  will 
nonetheless  be  slowed  down  (which  costs  valuable 
space,  but  recognizes  that  another  aircraft  might 
actually  be  just  ahead  of  the  simulated  flow). 
Inter-arrival  spacing  dips  whenever  an  aircraft 
that  happens  to  be  behind  the  TCSim  schedule  is 
followed  by  one  that  is  ahead,  because  the  ATSP 
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Fig.  4  Screen  snapshot  of  TCSim  simulating  merging  arrival  flows  along  the  charted  arrival  routes  shown  in 
Figure 


cannot  accelerate  the  lead  aircraft,  and  may  only 
slow  the  trail  aircraft  to  its  originally  scheduled 
arrival  time. 

•  Despite  differences  in  schedules,  if  TCSim  and 
the  algorithm  happen  to  agree  for  the  schedule 
of  two  successive  aircraft,  it  sometimes  starts  a 
new  sequence  of  aircraft  for  which  the  two  meth¬ 
ods  agree,  which  is  intuitive,  since  they  both  try 
to  assign  as  little  delay  as  possible,  once  the  con¬ 
straints  are  met. 

The  results  indicate  that  the  optimal  spacing  algo¬ 
rithm  can  provide  benefits  if  some  conditions  are  met. 
The  algorithm  would  then  produce  an  optimal  new 
schedule  for  an  arrival  flow.  First,  a  suitable  mecha¬ 
nism  for  generating  feasible  arrival  times  under  specific 
conditions  must  be  available.  This  is  an  instance  of  a 
class  of  problems  applicable  to  numerous  future  ATM 
concepts.  Second,  the  arrival  flow  must  be  reasonably 
well  established,  and  the  points  at  which  the  arrival 
times  are  generated  must  be  located  such  that  the 
times  become  available  to  the  algorithm  in  sequence. 


This  means  that  the  ATSP  can  immediately  issue  the 
advised  speed  to  the  aircraft.  Relaxing  this  require¬ 
ment  may  be  possible  by,  for  example,  repeating  the 
’speed  adjustment’  process.  Further  research  with  TC¬ 
Sim  will  explore  these  and  related  issues. 

Threshold  perturbation  intensity 

The  results  presented  above  also  illustrate  the  func¬ 
tionality  of  TCSim  for  which  we  want  to  use  this  al¬ 
gorithm:  for  a  given  perturbation  in  the  schedule  (due 
to  weather  or  other  factors),  which  can  be  simulated 
by  TCSim,  the  algorithm  enables  the  computation  of 
the  worst  case  scenario  under  the  optimal  action  of 
ATC.  In  future  implementations,  this  functionality  will 
thus  enable  us  to  compute  the  perturbation  threshold 
above  which  we  cannot  guarantee  a  given  spacing  re¬ 
quirement  at  the  destination  airport.  In  particular, 
we  are  interested  in  running  a  Monte  Carlo  simulation 
in  order  to  compute  the  threshold  covariance  of  the 
perturbation  above  which  airport  given  specifications 
become  impossible  to  meet. 
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Conclusion 


Fig.  5  Variation  of  the  maximum  minimum  time  sepa¬ 
ration  A  with  the  number  of  aircraft  which  have  crossed 
the  metering  fix.  This  value  is  obtained  by  averaging  20 
runs  and  will  be  obtained  by  a  Monte  Carlo  simulation 
in  the  final  version  of  the  paper. 


aircraft  number 


We  have  presented  some  of  the  simulation  capabil¬ 
ity  of  TCSim,  and  applied  it  to  the  Dallas/Ft.  Worth 
airport  airport.  As  a  building  block  of  our  method,  we 
have  derived  a  simple  dynamic  programming  based  al¬ 
gorithm  to  solve  for  the  maximum  minimum  spacing 
between  aircraft  upon  landing.  The  results  of  this  al¬ 
gorithm  were  compared  to  results  of  TCSim,  which 
are  not  only  based  on  spacing,  but  take  into  account 
other  factors  related  to  ATC.  We  have  shown  how 
these  results  will  be  used  in  future  to  characterize 
perturbation  thresholds  above  which  airport  spacing 
constraints  cannot  be  met  by  Air  Traffic  Control 
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Fig.  6  Example  of  spacing  computation  of  TCSim 
(solid  curve)  and  the  dynamic  programming  algorithm 
(dashed).  Each  point  of  the  curve  represent  the  A  com¬ 
puted  by  both  methods,  for  two  successive  aircraft. 
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